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ABSTRACT 

The stellar-dark halo mass relation of galaxies at different redshifts, M s (Mh,z), encloses relevant 
features concerning their physical processes and evolution. This sequence of relations, defined in the 
range < z < 4, together with average A cold dark matter (ACDM) halo mass aggregation histories 
(MAHs) are used here for inferring average M s growth histories, the Galaxian Hybrid Evolutionary 
Tracks (GHETs), where 'hybrid' remarks the combination of observational (M s ) and theoretical (Mh) 
ingredients. As a result of our approach, a unified picture of stellar and halo mass buildup, population 
migration, and downsizing of galaxies as a function of mass is presented. 

The inferred average M s growth histories (GHETs) of highest and lowest mass galaxies are defini- 
tively quite different from the average MAHs, A4(z), of the corresponding dark halos. Depending on 
how a given M^(z) compares with the mass at which the M s -to-Mh ratio curve peaks at the epoch z, 
Mh-p(z), two evolutionary phases are evidenced: (i) galaxies in an active regime of M s growth when 
Mh < Mh p , and (ii) galaxies in a quiescent or passive regime when Mh > Mh p . The typical M s at 
which galaxies transit from the active (star-forming) to the quiescent regime, M tran , increases with 
z, log(-Mt ran /-^©)~ 10.30 + 0.55z, making evident a population downsizing phenomenon. This result 
agrees with independent observational determinations based on the evolution of the galaxy stellar 
mass function decomposition into blue and red galaxy populations. The specific star formation rate, 
SSFR, predicted from the derivative of the GHET is consistent with direct measures of the SSFR 
for galaxies at different redshifts, though both sets of observational inferences are independent. The 
average GHETs of galaxies smaller than M tran at z = (M s » 10 10 3 M Q , M h « 10 118 M©) did not 
reach the quiescent regime, and for them, the lower the mass, the faster the later M s growth rate 
(downsizing in SSFR). The GHETs allow us to predict the transition rate in the number density of 
active to passive population; the predicted values agree with direct estimates of the growth rate in the 
number density for the (massive) red population up to z ~ 1. We show that ACDM-based models 
of disk galaxy evolution, including feedback-driven outflows, are able to reproduce the low-mass side 
of the M s -Mh relation at z ~ 0, but at higher z's strongly disagree with the GHETs: models do not 
reproduce the strong downsizing in SSFR and the high SSFR of low mass galaxies. 
Subject headings: cosmology: theory — galaxies: evolution — galaxies: haloes — galaxies: high- 
redshift — galaxies: star formation 
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1. INTRODUCTION two-point correlation function, galaxy group catalogs) to 

According to the popular hierarchical clustering A the theoretical Halo Mass Function (hereafter ffMF; for 

Cold Dark Matter (ACDM) scenario, galaxies form and ^nt reviews_ on all of these methods see Moster et al. 

grow inside evolving dark matter haloes. A central ques- 20 r 10 ; Behroozi et al. 2010; More et al. 2010, and more 

tion emerges then about what is the galaxy stellar mass, references m these papers). The latter approach is less 

M s , associated on average with a given halo of mass M h , dlrect than the f° mer * e > but ^ allows to cover larger 

i.e., the M s -M h relation. The change with redshift of mass ra ^ es ' and ^can be extended up to the redshifts 

this relation, M s (M h ,z), resumes the key astrophysical wh « e observed GS M Fs are reported (for recent results 

processes of galaxy stellar mass assembly in the context f e Conr °y & Wechsler 2009; Moster et al. 2010; Wang 

of the ACDM scenario & Jm § 2010; Be hroozi et al. 2010). 

With the advent of large galaxy surveys in the last ^ In Foster et *\ J 20 , 1 ?) B ehr00zi al ; ( 201 °; 

years, a big effort has been made in constraining the hereafter BCW10), the M s (M h ,z) functionality has been 

local M s -M h relation (1) directly by estimating halo constrained up to z « 4. In each case, different observa- 

masses with galaxy-galaxy weak lensing, with kinematics l 10nal data s ,f s im the GSMFs and different methods 

of satellite galaxies or with X-ray studies; and (2) indi- for s 1 tat if 1 tlca1 ^ ffSPmg halo masses to the galaxies were 

rectly by linking the observed statistical galaxy proper- used - The M s -M h relation at z ~ is similar m both 

ties (e.g., the galaxy stellar mass function GSMF, the wol] " s ( see also Baldr / et aL 2008 5 Dror y et aL ? 009 '> Gu0 

v 8 ' 8 J ' et al. 2010; Wang & Jing 2010). However, the change 

Electronic address: [firmani@merate.mi.astro.it,avila@astro.unam.mx| with Z of this relation, M s (M h ,z), particularly for Z>1, 

1 Also Instituto de Astronomia, Universidad Nacional Autonoma is different in both cases. 

de Mexico, A.P. 70-264, 04510, Mexico, D.F., Mexico Although the uncertainties in the inferences of 
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M s (Mh,z) are still significant, the current determina- 
tions can be used for preliminary explorations of the 
galaxy mass assembly process. An original approach has 
been introduced by Conroy & Wechsler (2009; hereafter 
CW09), who proposed a parametric form for GSMF as 
function of z constrained by both some observational re- 
ports of this function and the star formation rate (SFR)- 
M s relations to z ~ 1. Then, the cumulative GSMFs 
at each z were matched to the cumulative HMFs in or- 
der to infer the M s -Mh relations at different z's (abun- 
dance matching formalism, the simplest of the indirect 
methods; e.g., Marinoni & Hudson 2002, Kravtsov et al. 
2004, Vale & Ostriker 2004; see §2). Finally, the ob- 
tained M s -Mh relations at different z's were connected 
by using simple parameterizations of the average halo 
mass aggregation histories (MAHs) in order to infer the 
average stellar mass buildup of galaxies as a function of 
mass. 

The previous approach has the advantage (1) of be- 
ing flexible enough as to allow for explorations cover- 
ing large mass and redshift ranges, and (2) of providing 
a bridge between observational results and theoretical 
work. One of the bases of this approach is the statistical 
abundance matching formalism, which has shown to give 
robust results in agreement with direct inferences of the 
local M s -Mh relation or with other indirect inferences 
based on observations of the two-point correlation func- 
tions or galaxy group catalogs (see Moster et al. 2010; 
BCW10; More et al. 2010; Dutton et al. 2010a). 

From the purely empirical point of view, recent stud- 
ies are posing a new conception in our understanding of 
galaxy stellar mass assembly. For example, Drory & Al- 
varez (2008) used empirical GSMFs up to z ~ 5 (Drory 
et al. 2005) to infer the M s assembly of large galaxies, 
and in combination with available observations of the 
SFR-M S relation at different z's, they constrained the 
contribution of star formation (SF) and merging to stel- 
lar mass build up in galaxies. From this and other stud- 
ies (e.g., Bundy et al. 2006,2009; Hopkins et al. 2007; 
Pozzetti et al. 2010), the trend of downsizing (Cowie 
et al. 1996) is clearly shown in the sense that the mass 
at which the SFR starts to drop strongly decreases with 
time (from M s ~ 10 12 M at z > 4 to M s ~ 10 109 M Q 
at z ~ 0.5 according to Drory & Alvarez 2008). This re- 
sult combined with observationally inferred merger rates, 
led the latter authors to conclude that the (massive) 
red sequence is built up from top to down. This phe- 
nomenon was originally dubbed as 'archaeological down- 
sizing' (Thomas et al. 2005; see also for related results 
e.g., Daddi et. al 2004,2007; Drory et al. 2004; Bundy 
et al. 2005,2006; Conselice et al. 2007; Marchesini et 
al. 2009; Perez-Gonzalez et al. 2008 and more references 
therein), and it seems that its main driver are internal 
galaxy processes that efficiently quench SF in massive 
galaxies (e.g., Bundy et al. 2006; Peng et al. 2010). 

The latter result has been confirmed by recent decom- 
positions of the GSMF by galaxy type (based on color, 
SFR, morphology, etc.) up to z ~ 1 (e.g., Borch et al. 
2006; Bell et al. 2007; Drory et al. 2009; Ilbert et al. 
2010; Pozzetti et al. 2010; for a compilation of early ob- 
servations, see Hopkins et al. 2007). A systematic result 
obtained in these works is that blue/late- type (active) 
galaxies migrate to the red/early- type (passive) popula- 
tion involving a sequence of masses that decreases with 



time, a result in line with the archaeological downsizing 
phenomenon. The environment plays also an important 
role. Recent observational studies, where the GSMF has 
been divided not only by galaxy types but also by envi- 
ronment, show that the population migration happens 
more efficiently and earlier in denser environments (e.g., 
Peng et al. 2010). 

At lower masses the inferences of the M s assembly are 
difficultcd by the incompleteness limit of the samples as 
higher is z. However, at least for z<l, current studies 
show that the specific SFR (SSFR) of low-mass (blue) 
galaxies is surprisingly high, and higher than the SSFRs 
of more massive galaxies (e.g., Bauer et al. 2005; Noeske 
et al. 2007; Zheng et al. 2007; Bell et al 2007; Elbaz 
et al. 2007; Chen et al. 2009; Damen et al. 2009a,b; 
Santini et al. 2009; Oliver et al. 2010; Rodighiero et al. 
2010). This 'downsizing in SSFR' of low-mass galaxies 
seems to imply a delayed M s buildup as lower is the 
mass, something difficult to explain by the moment in 
ACDM-based models (for models, discussions, and more 
references see Noeske et al. 2007; Fontanot et al. 2009; 
Firmani, Avila-Reese & Rodrfguez-Puebla 2010). 

Based on the ACDM scenario and the new observa- 
tional inferences discussed above, a unified picture of 
stellar and halo mass buildup of galaxies of all sizes can 
be developed. By using M s -Mh relations constrained 
by BCW10 from z w to 4, we generate an analytical 
M s (Mh,z) relationship, which is continuous and differcn- 
tiable at any z. This sequence of M s -Mh relations in z 
is combined with average halo MAHs in order to: 

1 . calculate the average stellar mass buildup tracks of 
galaxies, M s (z), from z = 4 up to now (hereafter, 
galaxian hybrid evolutionary tracks, GHETs) along 
with their halo MAHs, M h (z); 

2. calculate the average SSFR histories of galaxies 
from the corresponding specific stellar mass growth 
rate histories, M s (z)/M s (z); 

3. determine the typical mass M tran at each z that 
marks the transition from active to passive galaxy 
population (in terms of M s growth) 

4. infer the transition rate in the number density of 
active to passive galaxies at each z. 

Our approach differs from CW09 in several aspects. 
These authors used observational inferences of the SSFR 
vs M s at different z's for constraining the evolution of 
their proposed GSMF (up to z ~ 1), which is a par- 
tial input of their approach for calculating M s (M\ 1 ,z). 
Instead, we use more recent and updated direct con- 
straints for this relationship (BCW10), which extends 
up to z w 4, and then predict the SSFR histories. In 
CW09 a simple parameterization for the average MAHs 
was used, while we generate them by using an extended 
Press-Schechter formalism that gives results similar to 
those of large N-body cosmological simulations. We dif- 
fer from CW09 also in several aims and results. 

The inferences of the quantities listed above constrain 
in general the multiple astrophysical processes partici- 
pating in the formation and evolution of galaxies. In par- 
ticular, we will compare the inferred evolutionary tracks 
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(GHETs) to predictions of standard ACDM-based semi- 
numerical models of disk galaxy evolution and explore 
whether the so-called downsizing in SSFR is really an 
issue for models. 

In §2, the method to calculate the GHETs is explained. 
The inferred average GHETs are presented in different 
diagrams in §3. In §§3.1 and 3.2 the evolution of the 
SSFR as a function of mass, the evolution of the transi- 
tion mass Mt ran , and the flow of active to passive galax- 
ies at each epoch, are shown and compared with direct 
observational estimates. The reliability of our approach 
and the results are discussed in §3.3. In §4, ACDM-based 
models of disk galaxy evolution are used to calculate evo- 
lutionary tracks for low-mass galaxies and see whether 
they are able to reproduce or not the GHETs, in par- 
ticular the SSFR downsizing phenomenon. A summary 
and our conclusions are given in §5. 

2. THE METHOD 

Our inference of the galaxy stellar mass buildup as a 
function of mass (called here GHETs) consists of two 
main steps: (1) we find an analytical functionality for 
M s (Mh,z)> and (2) after generating the halo MAHs, 
Mh(z), we combine them with M s {M^z) to infer the 
GHETs, M s (z). 

2.1. The Stellar Mass-Halo Mass Relations at 
Different Redshifts 

Our starting point is the determination of M s (Mh,z) 
up to the highest redshift as possible. As mentioned 
in Section 1, such a relationship has been recently con- 
strained up to z ~ 4 by using the statistical formalism of 
matching one-to-one the cumulative theoretical ACDM 
HMF to the observed cumulative GSMF at different 
z's: n h (> M h ) = n g {> M s ) (e.g., Moster et al. 2010; 
BCW10). Here we will use the results from BCW10, 
who started from an assumed parametric expression for 
Af s (Mh,z) and explored the parameter space, including 
uncertainties and sample variance, with a Monte Carlo 
Markov Chain technique to fit jointly the GSMFs ob- 
served at different z's (they actually carried out this 
analysis separately for two redshift ranges: 0<z< 1 and 
1 < z< 4). The authors have estimated and carefully taken 
into account all kind of uncertainties, which introduce 
significant scatter in the M s -Mh relations and affect even 
the shape of these relations (see §§3.3 for a discussion). 

The observational input used in BCW10 is the total 
local GSMF obtained in Li & White (2009) from Sloan 
Digital Sky Survey (SDSS) data, and the total GSMFs 
at higher redshifts (up to ~ 4) homogeneously obtained 
in Perez-Gonzalez et al. (2008). For the halo mass func- 
tion, BCW10 used the results from N-body cosmological 
simulations by Tinker et al. (2008). Since these HMFs 
refer only to distinct halos, a correction by the sub-halo 
abundance was introduced (see details in BCW10). The 
mass of the subhalo is fixed as that one at the time of ac- 
cretion into the halo, M acc , but the satellite galaxy's M s 
is associated with the current GSMF, i.e., it is implicitly 
assumed that the satellite's M s will continue to evolve in 
the same way as for centrals of halo mass M acc (the effect 
of fixing instead the satellite mass M s at the redshift of 
accretion has virtually no effect on the overall M s -Mh 
relation as checked in BCW10). Thus, the M s (M h ,z) 



relationship used here refers to the overall galaxy pop- 
ulation at each epoch, and this includes satellite galax- 
ies, though under the assumptions mentioned above for 
this sub-population. On the one hand, there are some 
hints that the M s -Mh relation of only satellite galaxies 
should not differ significantly from the overall relation 
(Wang et al. 2006). On the other hand, satellites are 
not a dominant population at any epoch, in particular 
at high redshifts. Therefore, even if its M s -Mh relation 
is quite peculiar, its effect over the average M s -Mh re- 
lation is expected to be negligible. For more details on 
the assumptions and role of uncertainties in the inference 
of the M s (Mh,2;) relationship we refer the reader to the 
original work by BCW10 (see also §§3.3). 

It is worth mentioning that BCW10 refer to their Ms- 
Mh relations as corresponding to "central" galaxies. This 
should be understood in the sense that for each (sub)halo 
is considered only the M s of one galaxy, the central one, 
but not in the sense that the satellite galaxy popula- 
tion has been excluded from the analysis (see above). A 
(sub)halo may contain more galaxies besides the central 
one (satellites). Therefore, the total stellar mass within 
the (sub)halo will be larger if we account for satellites. 
In BCW10 it was shown that the inclusion of satellites 
results in a "total" M s /Mh ratio sligthly larger than the 
"central" one at small masses but significantly larger at 
high (group and cluster) masses (BCW10; see also Yang, 
Mo & van den Bosch 2009). Some of this "excess" stel- 
lar mass is expected to end actually in the central galaxy 
due to satellite infall and/or tidal stripping of stars (other 
fractions of this stellar mass may end in the stellar halo 
or remain in the surviving satellites); this could modify 
the obtained overall M s -Mh relation 2 . In general, we do 
not expect that the overall M s -Mh relations inferred by 
BCW10 will be significantly affected by the (unknown) 
physics of satellite galaxies, at least for Mh 5 10 13 M at 
any epoch, and for any mass at z » 0. 

Figure Q] shows the M s -to-Mh ratio as a function of 
M h at z = 0.1, 1.0, 2.0, 3.0, and 4.0 reported in BCW10 
(solid, dashed, dot-dashed, dot-dot-dashed, and dot- 
dashed-dashed lines) with the respective 1<7 uncertainties 
(error bars, showed only for z = 0.1,2, and 4). Several 
features regarding galaxy mass assembly can be antici- 
pated from this figure: 

1) In the < z < 4 range the peak of the curves shifts 
from log(M hp /M Q )~ 11.8 to 12.5, while log(M s /M h ) p 
first declines from —1.5 at z « to —1.7 at z ~ 2, and 
then recovers the value —1.5 at z « 4. Taking into ac- 
count the uncertainty, the shape of the relations is con- 
served. Around log(Mh/M )^ 12. 5±0. 3 the curves cross 
each other. 

2) On the left side of the diagram, below the crossing 
masses (in particular from z = 1 to z = 0.1), for a given 
Mh, the M s -to-Mh ratio increases as z decreases, which 
reveals a strong late growth of M s for individual galaxies. 

2 According to the analysis of Yang et al. (2009), for halos 
smaller than < 10 13 1i -1 Mq at z ~ 0, the mass fraction of 
stars accreted by the central galaxy is negligible, indicating that 
the latter cannot have grown substantially due to the accretion of 
satellite galaxy's stars; rather their growth is dominated by indi- 
vidual evolution (in situ SF). For halos larger than Mb 10 13 
h _1 M0, the central galaxy M B may increase lately in a significant 
fraction by satellite accretion (dry mergers); another fraction of 
the accreted stars end in the stellar halo (intracluster stars). 
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Fig. 1. — Ms/Aih vs. Mji relations at z = 0.1, 1, 2, 3, and 4 from BCW10 with the respective uncertainties. The brown arrow shows the 
evolution of a passive galaxy (constant M s ) from z = 4, where its initial conditions coincide with the corresponding isochrone, to z = 0. 
The fact that such arrow reaches the 2 = isochrone shows that at these high masses galaxies evolve passively (see the text). The M B —M^ 
relations given by the M B (M\- l ,z) relationship proposed here are represented by dotted lines with the same color code as the BCW10 curves. 



Hence, this side is expected to be populated mainly by 
actively growing (blue) galaxies. 

3) Just the opposite happens on the right side, above 
the crossing masses, which means that for individual 
galaxies, M s brakes its growth while continues grow- 
ing. To illustrate this point, we plot in Fig. [I] the evo- 
lution (brown arrow) of a completely passive galaxy of 
constant mass M s « io 11 - 35 that starts in the isochrone 
z = 4 with Mh = 10 13 Mq and evolves until grows to 
» 1O 13 9 M at z = according to a typical ACDM MAH 
(see §§2.2) for this mass. If the galaxies of the sample 
that allow us to find the isochrones would be growing 
in M s , then the arrow would diverge from the z = 
isochrone. Therefore, the fact that the arrow reaches 
the z = isochrone exactly reveals that this side of the 
diagram is dominated by galaxies with passive evolution. 

The previous considerations confirm the idea that 
M s (M\ 1 ,z) resumes the key astrophysical processes of 
galaxy stellar mass assembly and therefore it contains in- 
formation about the evolution of individual galaxies. Our 
aim now is to explore the implications of the M s (Mh,z) 
relations shown in Fig. [T] As we are interested in the 
properties of M s (z) up to its second derivative (SSFR 
and the rate of change of SSFR), we cannot use the dou- 
ble parameterization in z given by BCW10 in the ranges 
< z < 1 and 1 < z < 4, respectively, because discon- 
tinuities arise around z = 1. Then, starting from the 
data of Fig. [TJ we introduce a modified parameterization 
that is continuous in z. The following equations, adapted 
from BCW10, summarize the model to be used here for 



< z < 4: 
Iog(M h (M.)) = 



The dependence on z is introduced in the parameters of 
equation (JTJ as: 

Iog(Mi(o))=Mi, + M M (a - 1), 
log(M Sj0 (a)) = M SA0 + M S)M (a - 1) + * (z) , 

P(a)=/3o+j3 a (a-l), (2) 
5(a) = S +5 a (a- 1), 
7(a) =7o + la (a - 1), 

where a = 1/(1+ z) is the scale factor. The function x( z ) 
controls the change with z of the peak value (ordinate) 
of the curves in Fig. [TJ If x ( z ) = 0, then the curve peak 
ordinate is independent of z. We chose x ( z ) m order 
to reproduce roughly the evolution of the peak found in 
BCW10. The first two Eqs. (2) control the position of 
the M s /Mh peak at each z (the value of Mh p [z]), while 
the last three ones control the shape of the M s / M^-M^ 
curves. We fix the parameters in Eqs. (1) and (2) ex- 
pecting (i) to reproduce well the BCW10 (M s /M h )-M h 
relations at the redshift extremes (z » and z ps 4, see 
Fig- UJi (h) to keep the generated relations at all z's 
within the la uncertainty given in BCW10, and (iii) to 
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TABLE 1 

Best parameters for the M- h (M a ,z) functionality 



Parameter 


BCW10 


This Work 




< z < 1 


< 2 < 4 




i/i ■)') — 

iu. <z_ 29 


10.70 


M s>0 , a 


U.OO_q 7 g 


—0.80 


A4i,o 


12.35lg;SS 


12.35 


M ha 


0.281°$ 


-0.80 


Pa 




0.44 


/3a 


o.i8ig-» 


0.00 


So 


u -°'-0.06 


0.48 


5a 


17+ - 42 


-0.15 


70 


1 56+ ' 12 


1.56 


7a 


n cm+0.15 
Z - 01 -1.S3 


0.00 



have the values of most of the parameters reasonably in- 
side the uncertainties of BCW10 (see their Table 2). The 
best set of parameters that we have found is reported in 
Table 1 (for comparison, we reproduce there also the best 
fit values from BCW10 for their < z < 1 case with free 
systematic parameters /i and re). For the function xi z )i 
we use: 

X (2) = -0.181z(l - 0.378z(l - 0.085z)). (3) 

The agreement of our model with the BCW10 results can 
be appreciated in Fig. [1] (dotted curves with the same 
color code as the results from BCW10 at z = 0, 1, 2, 3, 4). 

2.2. The dark halo mass aggregation history (MAR) 

The M s (Mh,z) relationship proposed in the previous 
section allows us to transform an average dark halo 
MAH, Mh(z), into an average stellar mass growth his- 
tory, M s (M h (z),z) (the GHET). Instead of using a pa- 
rameterization for the average MAHs as was done in 
CW09, here we calculate them by using the special Ex- 
tended Press-Schechter formalism developed in Avila- 
Recsc, Firmani & Hernandez (1998, see also Firmani & 
Avila- Reese 2000). Under the assumption that the pri- 
mordial density field is Gaussian, we calculate the overall 
mass distribution of progenitor halos at a given epoch 
Zi + i that will be contained in a halo of mass Mi at a 
later epoch Zj. Thus, after defining the cosmology and 
the mass power spectrum of fluctuations (as in BCW10, 
parameters very close to WMAP5 are used here, Ko- 
matsu et al. 2009), we start from a given mass at a given 
time (e.g., z — 0) and calculate its mass distribution of 
progenitors at a previous time step through Monte Carlo 
trials. The mass of the most massive halo (main pro- 
genitor) is used as the conditional for the next time step 
and so on. This way is obtained one realization of the 
MAH. For each M h given at z = 0, we calculate 20000 
random tracks (MAHs) and then we average among them 
to get the average MAH for this mass. In spite of the 
stochasticity, the MAHs present a systematical (hierar- 
chical) trend with mass: the smaller is M^(z = 0), the 
earlier is its mass assembly on average (Avila-Reese et 
al. 1998; van den Bosch 2002; Wechsler et al. 2002). 
Our average MAHs are not easy to describe by a simple 
parametrization, but they are in good agreement with 
those measured in the outcome of cosmological N-body 
simulations (e.g., Fakhouri, Ma, & Boylan-Kolchin 2010). 



3. GALAXY STELLAR MASS BUILDUP: AVERAGE 
EVOLUTIONARY TRACKS 

The GHETs in the M s vs. Mh diagram, calculated as 
explained in §2, are plotted with red thin solid lines in 
Figure [H The blue dashed curves show the M s -Mh re- 
lations for z — 0,1,2,3,4 from top to bottom (bottom 
to top) on the left (right) side. The GHETs in Figure 
[2] start from the z = 4 relation and end at the z = 
one. Two features of the GHETs in such a diagram are 
remarkable. (1) Below log(Mh/M Q )^ 11.6 the average 
slopes of the GHETs are d\ogM s /d\ogM\^ 4.5, greater 
than the slopes dlogM s /cflogM h ~ 2.3 of the M s -M h rela- 
tions in their low-mass side at any z. (2) The GHETs of 
massive halos (log(Mh/M Q ) > 12 at z = 0) stop to grow, 
happening this earlier, the more massive are the galax- 
ies. Such a behavior is due to the folding seen in the 
M s -Mh-;z surface that inverts the spatial distribution of 
the M s (Mh,z) curves going from the left to the right side 
of the diagram. The same properties may be derived 
from Figure |3j where the GHETs are shown as tracks 
that connect the M s / M^-M-^ relations. Here the folding 
is seen in the blue dashed lines as a maximum that shifts 
to the higher-mass side as z increases, making the lines 
cross each other. Both figures help to understand why 
the GHETs of massive halos attain a maximum. 

The GHETs ending at z = with masses 
8.4 <log(Ms/M Q ) < 11.6 are plotted with red solid lines 
in Figure |U The corresponding halo average MAHs are 
also plotted but, for comparative reasons, each one has 
been shifted vertically in such a way that each MAH coin- 
cides with its related GHET at z — (blue dashed lines). 
This figure shows in detail the difference between the av- 
erage stellar and the halo evolutionary track shapes. For 
M a (z=0)< IO 105 M , each GHET grows faster than the 
corresponding MAH since z = 4. For Af s (z=0)> IO 10 5 
Mq, as the system is more massive, the stellar assem- 
bly of the galaxy occurs earlier in time with respect to 
the corresponding halo. Our analysis shows conclusively 
that the stellar mass growth of galaxies deviates from 
the halo mass growth, specially for low- and high-mass 
systems. 

Figures [2] and |3] show that the more massive the galax- 
ies, the earlier their GHETs start to slow down their 
M s growth and at each epoch there is a typical halo or 
galaxy mass at which this growth stagnates; this typi- 
cal mass decreases with time. Such a situation can be 
interpreted as a transition at each epoch from active to 
passive galaxy populations, a transition at early epochs 
happens only for the largest galaxies but at later epochs 
involves gradually smaller galaxies (see §§3.2 for a more 
quantitative discussion). 

This transition from active to passive regimes is ob- 
served in the evolution of the blue and red components of 
the GSMF. In fact, observations show directly that, de- 
spite the stellar mass growth, the blue (late-type) galaxy 
component of the GSMF is not enriched with massive 
galaxies with time; on the contrary, within the uncer- 
tainty, it seems to become poorer in massive galaxies 
(e.g., Bundy et al. 2006; Bell et al. 2007; Pozzetti et 
al. 2010). Such a phenomenon is explained by a migra- 
tion of galaxies from the blue to the red component of 
the GSMF; the red (early-type) galaxy component of 
the GSMF is continuously fed with the migrated galax- 
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Fig. 2. — The M s -Mh relations (blue dashed lines) at z = 0, 1, 2, 3, and 4, from top to bottom (bottom to top) in the left (right) side, 
respectively, and the inferred average GHETs displayed from z = 4 to z = (red thin solid lines). The lower and upper dotted curves 
are model predictions for the evolution of two disk galaxies that end today with \o%,(M s /Mq) = 9.40 and 10.15, respectively (see §4 for 
details). From right to left, the crosses correspond to z = 0, 1, 2, 3, and 4. 



ies. Our result, in the light of such an interpretation, 
is remarkable because in our analysis the galaxy color 
bi- modality never has been taken into account. Studies 
where the GSMFs were decomposed into blue and red 
populations support our result in an independent way 
(e.g., Bell et al. 2007; Drory et al. 2009; Pozzetti et al. 
2010). 

Note that the prediction of a population transition and 
the typical mass that separates the two populations at 
each z depends on the position of the folding in the M s - 
Mh~z surface. This feature reveals a fundamental aspect 
of the GSMF evolution. Therefore, the accurate study 
of this feature deserves maximum efforts from the obser- 
vational point of view. 

The fact that the more massive are the systems, the 
earlier they finish their M s assembling, transiting from 
the active to the passive population, will be called here 
population downsizing (in the literature the term 'archae- 
ological downsizing' has been used to describe a related 
behavior of galaxies; see e.g., Fontanot et al. 2009, and 
the references therein). 

For systems less massive than M s « 10 105 M Q at z = 
(Mh « 10 12 M Q ), the active (growing) phase continues 
still at z = on average (the population downsizing has 
not happened), and the smaller is the system, the M s 
growth is more delayed and concentrated towards the 
present. The fact that the less massive is the system, 



the later happens its active phase of M s assembling, is 
called downsizing in SSFR (see e.g., Fontanot et al. 2009, 
Firmani et al. 2010 and the references therein). 

In the different diagrams plotted in Figures El |3l 
and SI we have introduced the corresponding evolution- 
ary tracks (crosses at z = 0,1,2,3, and 4 connected 
with dotted lines) for two (low-mass) galaxies calcu- 
lated by means of a ACDM-based model of galaxy evo- 
lution which will be explained in detail in §4. Their fi- 
nal stellar masses are log(M s /M Q ) = 9.4 and 10.1, re- 
spectively. We anticipate the rather different behavior 
of the model 'galaxian evolutionary tracks' (hereafter 
GETs) with respect to the GHETs. For example, the 
average slope of the former in the M s -Mh diagram is 
cflogM s /dlogMh=1.4, much less than the slopes of any 
GHET, and also less than the lower limit estimated above 
for the low-mass side of the M s -Mh relations at each z. 
This means that the model isochrones (defined by con- 
necting crosses of a given z) have a distribution and evo- 
lution opposite to that of the empirical Af s -Mh relations. 

3.1. Stellar mass growth rates as a function of mass 

The stellar mass buildup of galaxies may happen due 
to in situ SF or accretion of stars, mainly in dry merger 
events. Several theoretical (e.g., Mailer et al. 2006; 
Guo & White 2008), semi-empirical (Zheng, Coil Zc- 
havi 2007; CW09; Yang, Mo & van den Bosch 2009; 
Wang & Jing 2010), and observational (e.g., Bundy et 




Fig. 4. — Average GHETs (red solid lines) and the correspond- 
ing halo average MAHs (violet dot-dashed lines); for comparison 
purposes, each one of the latter ones was shifted vertically in such 
a way that M^(z=0) = M s (z=0). The lower and upper dotted 
curves are model predictions for the evolution of two disk galaxies 
that end today with log(M s /MQ) = 9.40 and 10.15, respectively 
(see §4 for details). 

al. 2006,2009; Bell et al. 2007; Drory & Alvarez 2008) 
pieces of evidence show that the former channel com- 
pletely dominates in low- and intermediate-mass galaxies 
at all epochs, while the latter may play a moderate role 



Fig. 5. — Evolution of the SSFR (see text for its definition) 
corresponding to the same GHETs shown in Figure [4] (solid red 
lines). From top to bottom, the curves with thick solid lines con- 
nect the tracks at the moment log(M s /Mo) equal to 9.5 (blue), 
10 (cyan), 10.5 (green), and 11.0 (red), respectively. The upper 
and lower black curves are model predictions for the evolution of 
two disk galaxies that end today with log(M s /M0) = 9.4 and 10.1, 
respectively (see §4 for details). 

for massive (M s ^10 11 Af Q , mainly red) galaxies at later 
epochs (z< 1). 
As a working hypothesis, we will assume that the stel- 
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lar mass buildup implied by our derived average GHETs 
is only due to in situ SF. Therefore, we calculate the 
SSFR as the time derivative of M s (z) divided by the cur- 
rent mass, M s (z)/M s (z), and divided by (1 — R), where 
R = 0.4 is the gas recycling factor due to stellar mass 
loss. The obtained SSFRs are then compared with di- 
rectly observed SSFRs as a function of M s and z. 

In Figure the evolution of the GHET-based SSFR is 
plotted (red solid line) for the same cases shown in Figure 
3J From top to bottom, the mass of the GHETs increases, 
except for the minor inversion (track crossing) present in 
the top right side of the diagram. For all the masses, 
SSFR decreases as z decreases. The overall behavior of 
decreasing SSFR with time agrees, by construction, with 
the known direct observational inferences of the cosmic 
stellar mass density history and its time derivative (e.g., 
Drory et al. 2005; Perez-Gonzalez et al. 2008), providing 
a comparison term to the observed cosmic SFR density 
history (Madau et al. 1996; see Hopkins & Beacom 2006 
for an extensive compilation of observational inferences). 

We test now if the GHET-based SSFRs histories as a 
function of mass are in agreement with the direct ob- 
servational measures. The observations of course do not 
refer to individual evolutionary tracks. What observers 
do is to determine the SSFRs of galaxy samples in dif- 
ferent redshift bins. A commonly reported result is the 
average of the SSFR in different z bins corresponding to 
a given range of M s (the same at all z's). In Figures [5] 
and [U from top to down, the thick (cyan, blue, green, 
and red) lines correspond to the SSFR inferred from the 
GHETs that in each z have log(M 8 /M©)=9.5, 10.0, 10.5, 
and 11.0, respectively. 

In Figure El we reproduce averages of SSFR obser- 
vational estimates at different z's for galaxies in four 
M s bins (centered in the same masses plotted for the 
GHETs). These estimates were taken for z ~ from 
Salim et al. (2007) excluding galaxies with AGN, and 
for higher z's (until z ~ 1), from Bell et al. (2007; only 
their Chandra Deep Field South estimates are plotted 
here). The data used to get these estimates refer to the 
population of blue (star-forming) galaxies (see Firmani 
et al. 2010 for more details on the choice of the obser- 
vational data). It should be stressed that the intrinsic 
spread around the plotted curves is rather large. 

For the largest mass bin (M s about 10 11 Mq, lower red 
curves), the GHETs imply clearly lower SSFRs than the 
directly measured SSFRs of star-forming (blue) galax- 
ies at z<l. This is expected since, as discussed above, 
the most massive galaxies show evidence to have early 
transited to the passive population. For the galaxy mass 
M s ~ 10 11 Mq, the transit to the passive, red sequence 
happens at z w 1.3 (see Fig. [5] below) as indicated by 
a vertical mark over the curve. Bell et al. (2007) re- 
ported the estimated SSFRs also for the red sequence 
galaxies. Their estimates for M s ~ 10 11 Mq (red circles) 
and 10 10 - 5 Mq (green circles) are plotted in Figure [51 
The directly measured SSFRs of massive (red) galaxies 
are very low, decreasing their values as z decreases in 
the same manner as our predictions show. While specu- 
lative by the moment due to the large uncertainties, the 
fact that the values of our inferred SSFRs are system- 
atically higher than the direct measures of red massive 
galaxies could be interpreted as an evidence of contribu- 
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Fig. 6. — Same curves connecting the GHETs at four values 
of Ms=constant (indicated inside the panel) as shown in Figure l5l 
(solid thick lines). The filled squares at z ~ are estimates from 
SDSS by Salim et al. (2007), and the small squares connected by 
dashed lines were inferred from Bell et al. (2007) for their sample 
of blue galaxies; in both cases the (average) masses are the same 
as the GHET-based masses. The filled circles are also inferences 
from Bell et al. but for their sample of red galaxies and only 
for log(M s /M )= 11.0 (red filled circles) and 10.5 (green filled 
circles). The vertical marks over the log(M s /Mo)= 11.0 and 10.5 
curves indicate the typical redshift at which these masses transit 
from active to passive phases. Note that the end of the 10 Mq 
GHET-based curve (when these masses are already in the passive 
phase) agrees well with the SSFR of red rather than blue galaxies; 
the end of the 10 10 5 Mq GHET-based curve is in between the 
red and blue galaxy samples. The dot— dashed line shows the curve 
l/[t[i(z) - 1 Gyr](l — R) corresponding to a constant SFR. 

tion of stellar accretion (dry mergers) to the stellar mass 
growth of these galaxies at late epochs (see Section 2.1 
for a discussion on this possibility). 

Recently, Damen et al. (2009a) presented estimates of 
the SSFR up to z w 2.8 for massive galaxies (M s > 10 105 
Mq) in the same Chandra Deep Field South used in Bell 
et al. (2007). As mentioned above, for redshifts z>l 
galaxies of M s » 10 11 Mq or smaller are on average in 
their active (SF-driven) growth regime. The results by 
Damen et al. are in rough agreement at these redshifts 
with our predictions. If any, the GHET-based SSFRs are 
slightly lower than the averages of the directly measured 
SSFRs. 

For the intermediate mass bin (M s k, 10 10 5 M©, green 
curves), the agreement between the average SSFRs in- 
ferred from the GHETs and those measured directly for 
the star-forming (blue) galaxies (Bell et al. 2007) is sat- 
isfactory, given the spread of data in SSFR. If any, the 
SSFR from the GHETs is slightly lower, in particular at 
low z's. Again, due to the population downsizing, ac- 
tive (blue) galaxies of masses Af s w 10 105 Mq transit to 
the passive (red) population at z « 0.35 (see Fig. [5]) 
as indicated by the mark over the green curve, so that 
their SSFRs at lower z's should be in between those mea- 
sured for the blue and red sequences. The agreement in 
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Fig. 7.— Same GHETs as shown in Figure [4] but in the SSFR 
vs M B (z) diagram (thin solid red lines). From bottom to top, the 
curves with thick solid lines connect the tracks at a given epoch 
(isochrones) : z =0 (cyan), 1 (blue), 2 (green), 3 (orange), and 4 
(red), respectively. The red crosses indicate the points where the 
GHET slope declines below -5. The lower and upper dotted curves 
are model predictions for the evolution of two disk galaxies that 
end today with log(M s /MQ) = 9.40 and 10.15, respectively (see §4 
for details). 

the mass range around 10 10 5 Mq with the Damen et 
al. (2009a) observational inferences up to z ~ 2.8 is also 
good, though the samples of these authors do not sepa- 
rate galaxies in early- and late-type ones. For z > 0.5 the 
slope of the SSFR(M S « 1O 1O ' 5 M )-(1 + z) relation for 
our prediction and for direct observations is about 3.3. 

For the lowest mass bins (M s ~ 10 10 and 10 9 5 Mq, 
blue and cyan curves), the agreement between GHET- 
predicted and directly-measured SSFRs is reasonable if 
takes into account (i) the large intrinsic spread of the ob- 
servational inferences (la width of ~ 0.5 dex around the 
mean in the relation log(SSFR)-logA/ s at z ~ 0, Salim et 
al. 2007; for higher redshifts, the typical reported intrin- 
sic ler widths are of ~ 0.2 — 0.3 dex, Noeske et al.2007; 
Zheng et al. 2007; Damen et al. 2009a,b); and (ii) the 
sample incompleteness, which increases as z is higher (for 
such small masses, this may bias the observational aver- 
age estimates toward larger values for z > 0.3 — 0.5). If 
any, the predicted SSFR(M s =const)-(l + z) curves for 
these low masses decrease slower (are shallower) than 
the curves inferred from direct observations for z < 1 
(see also Damen et al. 2009a). Better observational esti- 
mates for both the low-mass GSMF and SSFRs at high 
redshifts are needed to attain more conclusive results. 

The dot-dashed line in figure [6] shows the curve 
1/\Ph{z) - 1 Gyr](l — R) corresponding to a constant 
SFR. Galaxies above (below) this curve are currently 
forming stars at a rate higher (lower) than the past av- 
erage. 

In Figure the evolution of SSFR vs. M s for the 
same GHETs of previous figures is plotted (red solid 
lines). This plot remarks the fact that small galaxies 
keep high SSFRs and grow fast at late epochs, while very 
large galaxies had high SSFR and rapid growth in the 
remote past, being passive at later epochs. In Figure [3 



from bottom to top the thick (cyan, blue, green, orange, 
red) lines connect the individual tracks at the redshifts 
z = 0, 1,2,3, and 4 (isochrones), respectively. Each one 
of these isochrones corresponds to the average SSFR-Af s 
relations at a given z and can be compared with direct 
observational inferences of SSFR vs M s at different z 
bins. These inferences for samples at different epochs 
(from z ~ up to z ~ 2) show that galaxies (mainly 
star-forming ones) are ordered in the logSSFR-logM s di- 
agram along sequences of intrinsic widths ler< 0.3 — 0.5 
dex (e.g., Bauer et al. 2005; Salim et al. 2007; Schimi- 
novich et al. 2007; Noeske et al. 2007; Bell et al. 2007; 
Drory & Alvarez 2008; Damen et al. 2009a,b; Santini 
et al. 2009; Oliver et al. 2010; Rodighiero et al. 2010). 
The averages of these sequences are such that the smaller 
is M s , the higher the SSFR. However, due to the sam- 
ple incompleteness limit, which increases from a few 10 9 
Mq to several 10 10 Mq from z ~ to 1, respectively, 
these averages at low masses might be overestimated. 
The isochrones showed in Figure [7] are in general within 
the 1 a width of intrinsic spread of the observational 
SSFR-Af s relations. At low masses (where observational 
samples with measured SFRs are below the completeness 
limit), the isochrones tend to be below the direct obser- 
vational estimates of SSFR at the given epoch but at 
intermediate masses the agreement is rather good, i.e., 
at low masses, the GHET-based SSFR-M S relations tend 
to be shallower than those constructed from direct ob- 
servations. 

At the high mass end, our isochrones fall faster than 
the averages of the observed SSFR-M S sequences, which 
typically are reported for star-forming (blue) galaxies. 
Have in mind that our results refer to average trends 
and in this sense the average high mass galaxies in our 
case already transited to the passive (red) population. In 
fact, the SSFR-M S relation at a given z should be con- 
structed separately for active (blue) and passive (red) 
galaxies. Therefore, in the SSFR-M S diagram two se- 
quences would appear: at low masses a sequence domi- 
nated by blue galaxies and at high masses the sequence 
dominated by red galaxies. From the point of view of the 
GHET analysis, the characteristic mass that divides both 
sequences is related to the transition mass inferred from 
the GHETs at each z (see §§3.2 below). As seen in Fig. 7, 
the average slope of the passive sequence (high-mass side) 
at any z is steeper than that one of the active sequence 
(low-mass end). Interestingly enough, this is the trend 
actually revealed in those observational works where the 
SSFR-M S relations at different z's (out to z ~ 1.0 - 1.5) 
are plotted separately for red and blue galaxies (Bell et 
al. 2007; Oliver et al. 2010). 

In general, we conclude that the SSFRs at different 
z's and as a function of M s predicted on the basis of 
our GHETs are reasonably consistent with recent obser- 
vational studies based on direct measures of the SSFR 
of galaxies at different z's (especially those that divide 
galaxies into blue and red ones). This consistency is par- 
ticularly important because the observational informa- 
tion concerning SFR never enters in our scheme. There- 
fore, on the one hand, this result is an independent test 
for the whole approach presented here (it predicts the 
mean SSFRs of galaxies at different epochs consistent 
with direct observational estimates). On the other hand, 
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it supports our hypothesis that the stellar mass buildup 
given by the average GHETs is driven by in situ SF. Nei- 
ther does this mean that the growth of M s by accretion 
of stellar systems (dry mergers) may not happen, but 
it suggests that it is less dominant than the SF channel 
(our approach and the accuracy of current estimates of 
SSFR at different redshifts do not allow for quantitative 
inferences of the contributions of one or another channel 
of galaxy mass growth) . 

A natural question arises about the information the 
GHETs provide on the stellar initial mass function 
(IMF). We have outlined the agreement between our 
SSFR derived from the stellar mass growth and the di- 
rectly inferred SSFR. The latter inference is highly sensi- 
tive to the assumption of a universal IMF. In this sense, 
we can say that such an agreement is congruent with 
a universal IMF. Nevertheless, the uncertainties are yet 
large, making it premature to draw such a conclusion as 
sufficiently robust. 

Finally, the results mentioned above confirm and inte- 
grate into a unified picture two key facts of galaxy stellar 
mass assembly: (1) the downsizing in SSFR of low-mass 
galaxies and (2) the population downsizing of massive 
galaxies. The former will be discussed on the light of 
ACDM-based models in §4. The latter is quantified in 
the next section. 

3.2. Evolution of the characteristic mass that separates 
active from passive population 

Our method reveals the existence of a characteristic 
stellar mass at each z, M tT&n (z), at which the average 
GHET sharply declines its growth rate. We interpret this 
as a transition from the active to the quiescent regime 
of SFR, the onset of the migration from the star-forming 
blue to the passive red population (population downsiz- 
ing). We may estimate quantitatively M tran (z) by iden- 
tifying the mass at a given z corresponding to the GHET 
that started to strongly decrease its SSFR, for example, 
when the slope of the track in SSFR vs. M s (see crosses in 
Figure [7]) becomes steeper than —5 (the measured SSFR 
at this point increases from log(SSFR/Gyr _1 )= —0.8 at 
z ~ to ~ —0.5 at z ~ 3). The stellar Mtran calculated 
this way is plotted vs z in Figure |H] (diagonal crosses). 
The function: 

log(M tran /M Q ) = 10.30 + 0.55z (4) 
offers an approximate description of the results up to 
z ~ 2 (solid line). 

The focus now is to compare our result with inde- 
pendent observational pieces of evidence for a transition 
mass from active to quiescent/passive population as a 
function of z. Nowadays, as mentioned in the Introduc- 
tion, it is possible to follow the evolution of the early- 
(red) and late-type (blue) GSMFs separately. In one 
of the most recent and complete works, based on the 
zCOSMOS survey, Pozzetti et al. (2010) determined the 
mass where the early- and late-type GSMFs cross (dif- 
ferent estimators for these two populations are used, see 
the figure caption) from zwltozw0.2. Such a crossing 
mass, M cross , is interpreted namely as the typical mass of 
late-type galaxies migrating to early-type ones (see also 
Bell et al. 2007). 

The agreement between the Pozzetti et al. (2010) re- 
sults for M cross (circles with error bars connected by dif- 
ferent dashed lines) with our M tran (z) is satisfactory as 



can be seen in Figure [H Local estimates of M cross by 
Bell et al. (2003; filled triangle) and Baldry et al. (2004; 
filled square) are also plotted, as well as the law inferred 
by Drory & Alvarez (2008) for the mass above which the 
SFR as a function of M s begins to drop exponentially, 
M tran /M Q = 10 10 - 43 (l+z) 21 . Our result also agrees qual- 
itatively with other observational studies, with different 
definitions of the characteristic mass above which the 
passive (red) population of galaxies dominates in num- 
ber density (e.g., Bundy et al. 2006; Hopkins et al. 2007; 
Vergani et al. 2008). 

We may predict also the flow of galaxies transiting from 
active to passive populations at each z per unit of co- 
moving volume. At each redshift interval z, z + dz (time 
interval t, t+dt) there are GHETs transiting from the ac- 
tive to the passive regime, whose associated halo masses 
are in the interval Mh, Mh + t£M n . The transition rate in 
number density is then given by the abundance of halos in 
such a mass interval divided by dt, $(M n , z)xd\ogM\ 1 /dt. 
For the abundance of halos at each z we use the H M Fs 
given in Tinker et al. (2008) adapted to our cosmology 
and corrected to include sub-halos according to the func- 
tions given in BCW10. 

The result is shown in Fig. [S] The transition rate in 
comoving number density of active to passive galaxies up 
to z = 1 scatters around (1.0 — 5.5) x 10~ 4 gal Gyr -1 
Mpc -3 without any clear trend with z (the scatter is 
mainly due to the discretization of the MAHs in mass and 
z). At z ks 0.3, the rate is (3.7 ± 1.2) x 10" 4 gal Gyr" 1 
Mpc -3 . The obtained passive population growth rates 
can be compared with estimates reported in the same 
work by Pozzetti et al. (2010). They found an average 
growth rate in number density of the red population inte- 
grated above log(M s /M Q ) = 9.8 of 6.8(±1.2) x 10~ 4 gal 
Gyr -1 Mpc -3 for a redshift interval centered in z — 0.34 
(solid triangle with error bar in Fig. [SJ; this value is 
in reasonable agreement with our results. Note that the 
flow of galaxies in our case is related to a small mass 
range around the transition mass M tran given by the av- 
erage GHETs, while the direct observational determina- 
tion takes into account a large range of masses. The fact 
that the latter is only slightly larger than the former im- 
plies that indeed most of the galaxies becoming red are 
those of masses close to the average transition mass. 

3.3. Reliability of the results 

The results presented above are based mainly on the 
evolution of the M s -Mh relation. An important question 
is then how general and reliable is the M s (Mh,,z) relation- 
ship used here. A comprehensive analysis of the statisti- 
cal and systematical uncertainties in M a {M\ Xl z) and how 
they affect the shape of the M s -Mh relations at different 
z's has been carried out by BCW10. These authors sep- 
arate the uncertainty sources into three classes: uncer- 
tainties (i) in the observational inference of GSMF, (ii) 
in the dark matter HMF, and (iii) in the matching pro- 
cess arising primarily from the intrinsic scatter between 
M s and Mh. They found that by far the largest con- 
tributor to the error budget of the local M s -Mh relation 
comes from assumptions in converting galaxy luminosity 
into M s , which amounts to uncertainties of ~ 0.25 dex 
in the normalization of the relation (see the error bars 
in Figfl]). The contribution from all other sources of er- 
ror, including uncertainties in the cosmological model, is 
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Fig. 8. — Transition mass as a function of z (diagonal crosses), 
defined as the mass corresponding to the GHET that at a given 
2 started to brake its M B growth (the slope of the SSFR-M S rela- 
tion becomes steeper than —5). Masses above Mt ra n(z) are already 
in the passive phase (red quiescent population). The solid line is 
the dependence given in eq. The filled circles with error bars 
connected by lines are the crossing masses, M cross , at different z's 
determined in Pozzetti et al. (2010) as the cross of the GSMFs 
of the early— and late-type populations. The labels at the end of 
each curve indicate the different schemes adopted by the authors 
for the classification: MORPH1 and MORPH2 are morphologi- 
cal classifications using their ZH and MRS classification schemes, 
respectively; PHOT refers to a photometric classification scheme; 
SSFR1 corresponds to a classification into active and quiescent 
galaxies according to their SSFR. The dotted line is the law given 
in Drory & Alvarez (2008) for the mass above which the SFR as 
a function of mass begins to drop exponentially. The filled square 
and triangle are local determinations of the crossing mass of lo- 
cal blue and red GSMFs by Baldry et al. (2004) and Bell et al. 
(2003), respectively. 
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Fig. 9. — Transition rate in number density of active to passive 
galaxies at each z. The solid triangle with error bar is the esti- 
mate of growth rate in number density of red galaxies reported in 
Pozzetti et al. (2010) integrated above log(M s /M©) = 9.8. 



much smaller, ranging from 0.02 to 0.12 dex at z — and 
from 0.07 to 0.16 dex at z — 1. At high redshifts (z > 1), 
statistical uncertainties grow, becoming as significant as 
the systematic ones. On the other hand, the intrinsic 
scatter in M s given Mh and the random statistical error 
in M s inferences have a significant effect on the shape 



of the M s -Mh relation at the massive end (taken into 
account in BCW10). 

At z ~ 0, the results from most direct and indirect 
studies to infer M s (Mh, z ~ 0), in spite of the different 
methods and different local GSMFs used, agree well 
among them, especially if the uncertainties are taken 
into account (see BCW10 for a comparison). The sys- 
tematic shift of the (M s /Mh)--Mh curves to higher Mh 
as z increases (see Figs. [3]and[T]) is also a generic result 
obtained by all the authors (CW09; Wang & Jing 2010; 
Moster et al. 2010; BCW10). This implies that the 
trend reported here of downsizing in SSFR for galaxies 
less massive than M s « 3 1O 1O M0 is robust; the down- 
sizing can be respectively stronger or weaker depending 
on whether the low-mass end of the curves shifts more or 
less with z. As the completeness of the samples at low 
masses improves, the estimate of the GSMFs at their 
low-mass end will be better, and hence more accurately 
the evolution of the M s -Mh relation at lower masses will 
be constrained. 

The results of different authors may differ regarding 
the exact location of the M s -Mh (or M s / M^-M-^) curves 
at different z's, especially for M s > 10 10 5 Mq. For exam- 
ple, in Moster et al. (2010), the peaks in (M s /M h )- 
Mh significantly decrease with z in such a way that the 
curves intersect only at very high masses. This implies 
that the GHETs keep growing until z = even for large 
galaxies, the transition to passive population happening 
only for the most massive ones. In a similar way, CW09 
concluded from their analysis that out to z ~ 1 and for 
M s < 10 11 Mq there is not a characteristic mass at which 
the growth rate of galaxies is truncated. However, this 
is in conflict with independent observations that show 
instead strong population downsizing (see §§3.2 and Fig. 
i). 

In the case of BCW10, the peaks in (M s /M h )-M h 
slightly decrease out to z ~ 1 and for higher z's, in- 
crease attaining at z « 4 the same level as at z w 
(Figs. [T]and[3]). However, the uncertainties in the deter- 
mination of the (M s /Mh)-Mh curves are large and other 
behaviors with z are allowed. The relationship used here 
closely resembles the BCW10 results as discussed in §2. 

In order to probe the robustness of the conclusions pre- 
sented here, we have performed our analysis by using 
M s (Mh,z) relationships different to the one assumed in 
§2 but yet within the uncertainties of the BCW10 in- 
ferences. For example, we explored an M s (M} 1 ,z) rela- 
tionship such that the peaks of the (M s / M^-M^ curves 
remain at the same level for all redshifts (i.e., x(z)= in 
Eq. 2). The obtained GHETs are qualitatively similar 
to the ones presented here and Eq.Q slightly changes 
to log (M tran /M ) = 10.2 + 0.6z. It is noteworthy that 
the high mass GHETs in Fig. HI presenting a stagna- 
tion in their M s growth, show a shape a little sensitive 
to the parameter 8 a - By assuming S a — in Eq. (2), the 
most massive GHETs, after reaching a maximum stellar 
mass at approximately the same z as presented in Fig. 
0J decline in mass but not by more than 0.1 dex. 

We conclude that the main results reported here are ro- 
bust at least for the range of variations in the M s (Mh,z) 
relationship that remains within the current uncertain- 
ties in the inference of this relationship. 
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4. MODEL PREDICTIONS VERSUS HYBRID 
EVOLUTIONARY TRACKS 

Any model of galaxy formation and evolution within 
the hierarchical ACDM scenario should predict how does 
the M s -to-Mh ratio (the "galaxy formation efficiency" ) of 
individual galaxies evolve. Several complex astrophysical 
processes intervene in shaping the evolution of this ratio. 
Most of the semi-analytical and semi-numerical models 
-as well as the numerical simulations in that regards 
subgrid physics- introduce physical schemes for model- 
ing these processes and for reproducing the present- day 
M s -Mh relation (or the luminosity function) and other 
properties/correlations of galaxies. The question now is 
whether such models are able to reproduce the whole 
evolution outlined here. Different approaches to explain 
the population downsizing have been recently proposed, 
while the SSFR downsizing is just starting to be recog- 
nized as a problem (e.g., Fontanot et al. 2009; Firmani 
et al. 2010; Colin et al. 2010). 

In each one of the figures where the GHETs were pre- 
sented in this paper, we have plotted the GETs corre- 
sponding to two low-mass disk galaxy models that end 
today with log(M s /M Q ) = 9.40 and 10.15 (log(M h /M Q ) 
= 11.25 and 11.73; dotted lines). The models were calcu- 
lated with a semi-numerical method that self-consistently 
solves dynamical and hydrodynamical equations of: halo 
mass virialization given the MAH, disk formation in cen- 
trifugal equilibrium inside the growing halo, gravitational 
drag of the disk over the halo inner regions, SF triggered 
by gravitational instability and self-regulated by a verti- 
cal energy balance between energy input due to SNe and 
turbulent energy dissipation, and SN-driven mass out- 
flows (Firmani & Avila-Reese 2000; Firmani et al. 2010). 
The main assumptions implicit in the models are those 
of: spherical and cylindrical symmetries for the halo and 
disk, respectively, gas infalling on dynamical time scales 
('cold mode'), adiabatic invariance during halo contrac- 
tion, and detailed angular momentum conservation for 
the infalling gas. 

The baryon fraction is assumed initially as the univer- 
sal one, but the feedback-driven outflow reduces signifi- 
cantly this fraction. The parameters of the outflow model 
are fixed in such a way that namely the (low-mass) Ms- 
Mh relation at z ~ is reproduced. In Firmani et al. 
(2010) we have shown that this happens only in the case 
of the so-called energy-driven outflow and for a high SN 
energy-transfer efficiency 3 . Because the exploration in 
this paper is only at the level of average trends, the halo 
MAHs used in the models are the averages among 20000 
realizations for each given mass, and the used halo spin 
parameter A was fixed to 0.03 at all z's, a value slightly 
smaller than the mean of z — relaxed halos in N-body 
cosmological simulations (e.g., Bett et al. 2007). On 
the one hand, the N-body simulations show that A does 
not significantly changes with time for halos growing in 
the accretion mode (e.g., Peirani, Mohayaee & de Freitas 
Pacheco 2004; D'Onghia & Navarro 2007). On the other 
hand, our models show that the radius evolution of disk 
galaxies agrees with observational inferences when the 

3 If the possibility of late re-accretion of the lost gas is taken 
into account, then in Firmani et al. (2010) it was found that the 
momentum— driven outflow is the one that best can help to repro- 
duce the 2 = M B -Mh relation. 



gas A parameter is roughly constant in time (Firmani & 
Avila-Reese 2009). 

The GETs in Figure 0] show that the stellar mass 
growth is roughly proportional to the halo MAH as z 
decreases. The model GETs and the GHETs in the 
M s -Mh diagram (Fig. [2]) are very different; the average 
slopes of the GETs are cflogM s /cflogMh~1.4, much less 
than the slopes of the corresponding GHETs. The model 
M s -to-Mh ratio in Figure [3] decreases toward the past 
due mainly to the feedback-driven outflow effect, but 
slightly, while the GHETs of the same masses evidence 
a very fast decreasing. The M s -M h and (M s /M h )-M h 
relations (Figs. [2] and El respectively) predicted by the 
model GETs at low masses change with z (M s slightly 
increases with z for a fixed M^) in an opposite direc- 
tion with respect to the empirically inferred observations. 
This discrepancy represents a serious problem for all the 
theory behind the models. 

The dramatic differences between the GETs and 
GHETs are also seen in the SSFR evolution plots (Figs. 
[5] and [7]). We remark here the downsizing aspect. The 
GET-based SSFRs of the two galaxies evolve almost par- 
allel (with a law SSFRcx (l+z) 2 - 2 ~ 2 - 3 ), being the SSFR 
of the larger galaxy slightly higher since high z's (weak 
upsizing). Instead, the GHET-based SSFRs cross each 
other at high z (see the upper side of Figure [5J, in such a 
way that larger galaxies have higher SSFR at very high 
z's, but at lower z's, the situation is inverted in such 
a way that as z — > 0, the smaller is the galaxy, the 
higher the SSFR (downsizing). The comparison of the- 
oretical models with the GHET-based SSFRs suggests 
that some astrophysical mechanism not considered up to 
now should systematically delay the onset of SF activity 
the smaller the halo is. 

The low SSFRs (early stellar mass assembly) of small 
galaxies predicted by the models, which include energy- 
driven outflows, are actually generic to all ACDM-based 
models. This issue has been reported in different ways 
by the semi-analytic models. Most of these models 
show that the stellar population of small galaxies (M s ss 
10 9 — 10 10 ' 5 M Q ) is assembled too early, becoming these 
galaxies older, redder, and with lower SSFRs at later 
epochs than the observed galaxies in the same mass range 
(e.g., Somerville et al. 2008; Fontanot et al. 2009; Santini 
et al. 2009). The problem is both at the level of satel- 
lite and central galaxies. For example, it was shown that 
models typically overpredict the observed stellar popula- 
tion ages (Pasquali et al. 2009) and stellar masses (Liu et 
al. 2010) of central galaxies in low mass halos. By means 
of a galaxy evolutionary model similar to ours, Dutton, 
van den Bosch & Dekel (2010b) also find that the SSFR 
of model low-mass galaxies is below the average of ob- 
servations, specially at redshifts z ~ 1 — 2, though these 
authors conclude that their models are able to reproduce 
roughly the main features of the observed SFR sequence. 

Finally, it should be said that N-body + hydrodynam- 
ics simulations of individual low mass galaxies in a cos- 
mological context also face the issue of too low SSFRs (as 
well as too high stellar and baryonic mass fractions; e.g., 
Colin et al. 2010). In Firmani et al. (2010) we explored 
the inclusion in our semi-numeric models of re-accretion 
of the ejected gas for a broad range of possibilities. As 
expected, the SSFR of galaxies increases in general due 
to re-accretion but this increasing goes in the opposite 
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direction of downsizing: only moderately for low-mass 
galaxies and too much for the larger galaxies. 

5. SUMMARY AND CONCLUSIONS 

Empirically inferred M s -M n relations from z ~ to 
z ~ 4 were connected with average ACDM-halo MAHs 
in order to infer the corresponding individual average 
stellar mass growth histories, called here GHETs. We 
have adopted an M s (M n ,z) relationship continuous in the 
— 4 redshift range and in agreement with the inferences 
by BCW10, who used the technique of jointly matching 
the abundances of observed GSMFs to the theoretical 
HMFs at different z's. The main results obtained here, 
which allow to establish a unified description of average 
galaxy evolutionary tracks for a large range of masses, 
are as follows. 

• The folding in the Mg-M-^-z surface, reflected mainly 
as peaks in the (M s /M n )-Mh diagram that shift to 
higher masses as z increases, introduces an important 
feature in the behavior of the GHETs. For masses much 
smaller than the peak mass at a given z, Mh p (z), the 
GHETs are in an active growth phase, whereas masses 
close to or greater than M np (z) are in a quiescent or 
completely passive (stagnated) phase (Figs. [2] and [3]). 
Therefore, at each z there is a characteristic stellar mass 
at which on average galaxies slow down their growth 
and transit from an active (star-forming blue) popula- 
tion to a passive (red) population. This transition mass 
decreases with time, giving rise to a phenomenon called 
here 'population downsizing'. 

• Galaxies less massive than M np at z = (M s < 10 105 
M Q , M h <10 120 M Q ) have still growing GHETs. Be- 
sides, the lower the mass, the faster the later M s growth, 
due likely to a delayed and lately active SF phase, a phe- 
nomenon called 'downsizing in SSFR '. 

• The shapes of the average stellar and halo mass as- 
sembling histories are quite different (Figure |4j. For 
galaxies that at z = have M s < 10 10 5 M , their dark 
MAHs at later epochs grow slightly slower as smaller is 
the mass, while their stellar GHETs grow much faster. 
For M s (z = 0) > 10 10 - 5 M Q , the larger the galaxy, the 
earlier its GHETs attain a stellar mass stagnation; until 
this epoch, the more massive the galaxy, the faster the 
M s growth with respect to the corresponding Mh growth. 

• By neglecting any stellar mass growth by accretion 
(dry mergers) , a GHET allows to find the corresponding 
SSFR history, SSFR(z), by calculating M s /M s and cor- 
recting by the stellar mas-loss recycling factor R. The 
inferred SSFRs at different z's and as a function of M s 
(Figs. [6j and [7|) are reasonably consistent with re- 
cent observational studies based on direct measures of 
the SSFR of star-forming galaxies at different z's. The 
GHET-based SSFRs corresponding to masses that at a 
given z are larger than the transition mass M% Tml (z), are 
much lower than the measured SSFRs of (rare) luminous 
blue galaxies, but agree or are even slightly larger than 
the measured SSFRs of red galaxies. The overall con- 
sistency between the GHET-based SSFR-Af s relations at 
different z's with direct inferences of these relations sug- 
gests that the accretion of stellar systems (mergers) plays 
a minor role in the M s assembling of galaxies, excepting 
perhaps those that transited to the passive sequence (the 
most massive ones). 



• By using the SSFR vs. M s evolutionary tracks, we 
calculated the characteristic transition mass at each z, 
M tran (z), above which the average GHET starts to sig- 
nificantly decrease its growth rate, transiting from the 
active to the passive galaxy population. The result 
is roughly described by the relation log(M tr an/M Q ) = 
10.30 + 0.55z (at least up to z ~ 2). This result agrees 
with recent observational determinations of the evolution 
of the mass, at which the early- and late- type GSMF 
components cross each other, from z w 1 to the present 
(Fig. U). 

• We determined also the transition rate in number 
density of active (blue) to quiescent (red) galaxy popula- 
tion. At z w such a rate is i()- 3 - 4±0 - 2 gal Gyr -1 Mpc~ 3 
and up to z = 1 the rates are within (1.0 — 5.5) 10~ 4 gal 
Gyr -1 Mpc~ 3 , in good agreement with direct observa- 
tional inferences of growth rate in number density of red 
galaxies. 

• We further explored whether ACDM-based models 
of galaxy evolution are able to predict galaxy evolution- 
ary tracks (called here 'GETs') in agreement with the 
GHETs, in particular at low masses where the downsiz- 
ing in SSFR is evidenced by the GHETs. We have shown 
that while the models with SN energy-driven outflows 
are able to reproduce the local M s -M n relation (see also 
Firmani et al. 2010), they fail in reproducing this relation 
at higher z's. The difference between GETs and GHETs 
at low masses is rather large: the former ones show a fast 
decrease in SSFR with time almost independent of mass, 
while for the latter ones, the lower the mass, the slower 
the late decrease of SSFR (downsizing in SSFR). 

From these results, we conclude that the general de- 
scription of stellar mass buildup of galaxies provided by 
our average GHETs appears rather successful: the pre- 
dicted SSFR histories as a function of M s and the pre- 
dicted transition mass M tran as a function of z, as well 
as the transition flux at different z's, are roughly con- 
sistent with direct (but yet limited) observational esti- 
mates of all these quantities. Furthermore, our analysis 
reveals nicely the existence of a galaxy bi-modality re- 
lated to the M s growth activity: as z decreases, smaller 
and smaller masses transit on average from the active 
(blue star-forming) to the passive (red) population. 

An important ingredient of our approach is the connec- 
tion of ACDM halos to the observed galaxies. Therefore, 
the underlying hierarchical ACDM scenario is successful 
in the sense that the predictions of the model are in agree- 
ment with independent observations. The average stellar 
and dark mass buildup of galaxies were found to be signif- 
icantly different, specially for low and high masses. This 
result is in qualitative agreement with CW09 in spite of 
the differences in the data, method, and redshift range 
between their and our work. However, the results from 
both works are different at a quantitative level, which is 
seen, for example, in the different conclusions regarding 
the existence of a characteristic mass at each epoch at 
which the SSFR of galaxies is truncated. Note that the 
SSFR and the determination of this characteristic mass 
(Mtran) are related to the first and second derivatives of 
the M s growth, respectively. 

Our results put in a unified picture both the popula- 
tion downsizing (related to galaxies with M s > 3 1O 1O M 
at z = 0) and the downsizing in SSFR (related to smaller 
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galaxies). The challenges are now to explain the pro- 
cesses that produce: (1) the transition from active to 
passive regimes in M s growth associated with a sharp ces- 
sation of the SFR in massive galaxies, (2) the fact that 
the typical mass of SSFR truncation, M tTan , decreases 
with cosmic time (population downsizing), and (3) the 
fact that the smaller is the halo, the more is delayed the 
galaxy's M s growth (downsizing in SSFR). 

Items (1) and (2) seem to find partial explanations in 
aspects related to the same dark halo clustering (envi- 
ronment) evolution (e.g., Neistein et al. 2006) as well 
as to the introduction of feedback processes due to the 
AGNs of massive galaxies (e.g., Silk & Rees 1998; Kauff- 
mann & Haehnelt 2000; Granato et al. 2004; Cattaneo 
et al. 2005; de Lucia et al. 2005; Bower et al. 2006; 



Croton et al. 2006). Item (3), as shown here, is a sharp 
and less understood problem. Later re-accretion of the 
gas ejected by galaxies due to SN-driven outflows does 
not account for a solution to this issue (Firmani et al. 
2010). A better understanding of the astrophysical pro- 
cesses intervening in galaxy formation and evolution as 
well as the role of environmental processes is necessary. 
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